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FERMI LAT OBSERVATIONS OF SUPERNOVA REMNANTS INTERACTING WITH MOLECULAR CLOUDS 
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ABSTRACT 

We report the detection of 7-ray emission coincident with four supernova remnants (SNRs) using data from 
the Large Area Telescope on board the Fermi Gamma- ray Space Telescope. G349.7+0.2, CTB 37A, 3C 391, 
and G8. 7-0.1 are SNRs known to be interacting with molecular clouds, as evidenced by observations of hy- 
droxyl (OH) maser emission at 1720 MHz in their directions. SNR shocks are expected to be sites of cosmic 
ray acceleration, and clouds of dense material can provide effective targets for production of 7-rays from tt"- 
decay. The observations reveal unresolved sources in the direction of G349.7+0.2, CTB 37 A, and 3C 391, and 
a possibly extended source coincident with G8. 7-0.1, all with significance levels greater than lOcr. 
Subject headings: gamma rays: ISM — ISM; individual (G349.7+0.2, CTB 37A, 3C 391, G8.7-0.1) — ISM: 
supernova remnants 



1. INTRODUCTION 

The case for cosmic ray acceleration at supernova rem- 
nant (SNR) shocks has found increasing support in observa- 
tional evidence, both from the thermal and non-thermal com- 
ponents of the emission from these objects. The production 
of relativistic electrons has been clearly established from the 
non-thermal X-ray em ission detected from several shell-typ e 
SNRs, e.g., SN 10 06 (iKovama et al.1 



1998), 



RX J1713.7 -3946 (IKovama etal.1 II 997L ISlaneetar.l Il999h. 
and Vela Jr. ( lAschenb ach"1998l: ISlane et al.ll2001l) . Addition 
ally, the dynamical properties of some SNRs suggest that a 
significant fraction of the explosion energy is going into the 
acceleration of i ons to relativistic sp eeds, like in the case of 
IE 0102.2-7219 (iHughes et al.l2000l) . X-ray measurements of 
the relative positions of the forward and reverse shocks, and 
contact discontinuity in Tycho's SNR also indicate modifica- 
tion of t he remnant dynami cs by the acceleration of cosmic- 
ray ions (IWarren et al.l2005b . 7-ray emission in the TeV range 
from some SNRs, observed with ground-based telescopes, 
seems to confirm that particles are indeed being accelerated 
to energies approaching the knee of the cosmic ray spectrum 
jEnomoto et al."2002'; 'Aharo nian et al.l 12007 at iKatagiri et all 
12005; Aharonian et al..,2007bi) . However, there is growing de- 
bate as to whether the origin of this high energy emission 
is hadronic or leptonic in nature, as reviewed in iReynoldsl 
(I2OO8). 

Pion production, followed by tt"— >77 decay, is a natural 
result of the interaction of relativistic ions, like those gen- 
erated through diffus ive shock accelera t ion at SNR shocks , 
with dense material dPrurvet al.l 119941: iBaring et all 119991) . 
Hence, SNRs known to be interacting with molecular clouds 
present particularly interesting targets for the dete ction of 7- 
ray em issi on of hadronic origin, as suggested by iFrail et aP 
(119961) and lClaussen etal.1 ( 119971) . The detection of hydroxy 1 
(OH) masers at 1720 MHz is a very useful tool for finding this 
type of interactions, as they indicate t he presence of shocked 
clouds of dense rnolecu lar material (IWardle & Yusef-Zadehl 
IM)llHew^tterd1l2009l) . 

Observations in the MeV to GeV range are especially im- 
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portant in order to identify signatures of 7r°-decay, because 
the photon spectrum from this pr ocess has a dist ribution that 
peaks approximately at 70 MeV ( lRevnoldsll2008l) . The Ener- 
getic Gamma-Ray Experiment Telescope (EGRET), on board 
the Compton Gamma-Ray Observatory, found a host of dis- 
crete 7-ray sour ces in the GeV ban d possibly associated with 
Galactic SNRs ( iTorres et al.l 12003*). However, due to poor 
spatial resolution and other instrumental constraints, it has 
proven impossible to determine unambiguously the origins of 
these EGRET sources. The Large Area Telescope (LAT), on- 
board the Fermi Gamma-ray Telescope has opened a new win- 
dow for studying SNRs in the GeV energy range, due to its 
improved sensitivity and resolution. For example, the detec- 
tion of an extended LAT sourc e coincident with SNR W51C 
(G49. 2-0.7), was reported by ( lAbdo et alj|2009d) . and their 
study suggests that 7r°-decay emission is the dominant com- 
ponent of the 7-ray flux observed in its direction. The obser- 
vation of OH maser (1 720 MHz) emission in the direction of 
this SNR, reported by iGreen et al] ( 119971) . suggests that it is 
interacting with a cloud of dense material. 

SNRs G349.7H-0.2, CTB 37A, 3C 391 and G8.7-0.1, are 
known to be interacting with molecular clo uds, as evidenced 
by OH masers observed in th eir directions (IFrail et al.l 1 19961 : 
iHewitt & Yusef-Zadehll2009l) . In this paper we report the de- 
tection of four LAT sources coincident with these SNRs, and 
on the results of our detailed study of 7-ray emission in their 
regions. After a description of how the LAT data are analyzed 
in this work in Section 2, we present the results of the obser- 
vations and put them in context for each individual object in 
Section 3. A discussion of the results is given in Section 4. 

2. OBSERVATIONS AND DATA TREATMENT 

The LAT 7-ray detector is the main instrument on board the 
Fermi Gamma-ray Space Telescope, which was launched on 
2008 June 1 1 . The LAT is a pair-conversion, imaging, large 
field-of-view (FoV) 7-ray telescope, that covers a wide en- 
ergy range, from 20 MeV and up to 300 GeV. High precision 
tracking of electrons and positrons resulting from pair con- 
version allow the LAT to reconstruct the direction of 7-rays 
incident on the detector, and measurements of the energy of 
the subsequent electromagnetic showers are performed in its 
calorimeter The FoV is ~ 2.3 sr and the point-spread func- 
tion (PSF) varies largely with photon energy, and improves at 
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higher energies. The converter-tracker device has 16 planes 
of high-Z material, which are divided in two regions. The 
front section contains the first 12 layers which are thinner and 
optimize the PSF, and the back region has 4 thicker planes, 
which maximize the effective area at the expense of a lower 
angular resolution relative to the front section. Further de- 
ta ils about the instrum ent and data calibration are described 
in lAtwood et"an (l2009l) . 

The LAT, while operating on "sky survey" mode, observes 
the entire sky every two orbits (approximately 3 hr). In this 
work, 15 months of LAT data (from 2008 August until 2009 
November) are analyzed. Only events belonging to the dif- 
fuse class, which r educes the residual ba ckground rate as ex- 
plained in detail in I At wood et al.l ( 120091) . have been selected 
for this study. The updated instrument response functions 
(IRFs) used are called "Pass6 version 3", which were devel- 
oped using in-flight data and consider pile-up and coincidence 
effect s not included in the pre-launch analysis (iRando et alJ 
120091) . The systematic uncertainties of the effective area, for 
the IRF used, are energy dependent; 10% at 100 MeV, de- 
creasing to 5% at 5 60 MeV, and increasing to 20% at 10 GeV 
(lAbdo et al.l 120091) . Additionally, only events coming from 
zenith angles smaller than 105° are selected in ord er to reduce 
the co ntribution from terrestrial albedo 7-rays (lAbdo et alJ 
120093) . 

The 7-ray data in the direction of G349.7H-0.2, CTB 37A, 
3C 391 and G8. 7-0.1 are analyzed using the Fermi Science 
Toolfl In order to study spatial distribution, position and 
spectral characteristics, the maximum likelihood fitting tech- 
nique is employed using gtlike. The likelihood statistical ap- 
proach is used on LAT data because of the low detection rates 
and the extent of the PSF, and it estimates the best model 
parameters by maximizing the joint probability for the data 
given the emission model (iMattox et al.l 119961) . The emis- 
sion models used in gtlike place individual sources, with par- 
ticular spectral properties, at fixed positions, and include a 
Galactic diffuse component resulting from the interactions 
of cosmic rays with interstellar material and photons, and 
an isotropic one that accounts for the extragalactic diffuse 
and residual backgrounds. In this work, the mapcube file 
gll_iem_v02.fit is used to describe the 7-ray emission from 
the Milky Way, and the isotropic component is modeled using 
the isotropic_iem_v02.txt table. 

2.1. Spatial Analysis 

Spatial analysis is performed using 7-ray data in the energy 
range 2-200 GeV, and converted in t\\t front section, in order 
to improve the angular resolution. The 68% containment ra- 
dius angle for normal incidence/ronf-selected photons in this 
energy band is < 0°.3. To analyze the neighborhoods of the 
SNRs studied, count maps of the 10° x 10° region centered 
on the region of interest (ROl) are constructed, and smoothed 
by a Gaussian kernel to match the PSF of the events selected. 

In order to determine the detection significance and the pre- 
cise position of the sources associated with the objects of in- 
terest. Galactic and isotropic backgrounds are modeled and 
test statistic maps are constructed using gttsmap. The test 
statistic is the logarithmic ratio of the likelihood of a point 
source being at a given position in a grid, to the likelihood of 
the model without the additional source, 21og(Lps/Lnuii). The 
resolution of these maps, set by the size of the grid used, is 

^ The Science Tools package v9rl5p2, and related documentation, is dis- 
tributed by the Fermi Science Support Center at http://fermi.gsfc. nasa.gov /ss^ 



0°.05. 

2.2. Spectroscopy 

The study of the spectral energy distribution (SED) charac- 
teristics of each of the sources associated with G349.7-I-0.2, 
CTB 37A, 3C 391, and G8.7-0.1, is performed through like- 
lihood analysis with gtlike. Events converted in both/ro«r and 
back sections, and in the energy range 0.2-51.2 GeV, are in- 
cluded in the study. The lower energy bound is selected to 
avoid source confusion, to avoid the rapidly changing effec- 
tive area of the instrument at low energies, and because of the 
large uncertainty below 0.2 GeV related to the Galactic dif- 
fuse model used. 

The spectral analysis is performed using gtlike to model the 
flux at each energy bin and estimate, through the maximum 
likelihood technique, the best-fit parameters. Background 
sources fro m the 3 month Fe rmi LAT Bright Gamma-ray 
Source List dAbdo et alj2009al) . and dete cted pulsars from th e 
1 1 month Fermi LAT Pulsar Catalogue ( lAbdo et alj|2009bl) . 
within a 10° radius of the sources of interest are modeled. 
Sources detected in the test statistic map of each field with 
significance greater than 5(t are also included in the back- 
ground model. Additionally, the analysis was repeated us- 
ing all background sources in the recently released one-year 
Fermi LAT First Source Catalog (IFGL) in each regiorQ, with 
no significant effect on the results. The normalization of the 
Galactic diffuse component is left free at each energy step. 
As an addition to the statistical uncertainties associated with 
the likelihood approach, and the systematic errors related to 
the IRF as mentioned above, the uncertainty of the under- 
lying Galactic diffuse level is considered. The treatment of 
this error in the evaluation of the systematic uncertainties is 
performed by artificially changing the normalization of the 
Galactic backgroun d by ±3%, similarly to the analysis in 
lAbdo et all (120093) . 

3. RESULTS AND ANALYSIS 

3.1. G349.7+0.2 

G349-I-0.2 has one of the highest radio and X-ray surface 
brightnesses in the Milky Way. Its radio image is character- 
ized by a small diameter (approximately 2'. 5) partial shell, 
somewhat enhance d toward the southern and eastern limbs 
dShaver et al.lll985"l) . Five OH maser emission spots (at 1720 
MHz) were detected in the direction of this SNR, which 
suggests that it i s interacting with a dense molecular cloud 
dFrail et al.l 1 19961) . The velocity of these maser lines places 
the regions of shocked material at a distance d ~ 22.4 kpc. 
The characteristics of the emission detected from molecular 
transitions of '^CO, '^CO, CS, HCO+, H2CO, HCN and SO, 
at similar velocities to those of the OH masers , are further ev- 
idence of molecular cloud-shock interaction dLazendic et alj 
l200llRevnoso & Mangumll2001l) . 

In the X-ray band, the morphology of the remnant ob- 
served with ASCA and Chandra is similar to that from the 
radio, and the high photon absorption of its X-ray spectrum 
(A^H « 7 X 10^^ cm"^) is consistent with the distance esti- 
mate from maser emission dYamauchi et alJll998l;ISlane et alj 
l2002t lEazendic et al.l 120051) . In the analysis of the Chandra 

The data for the 1451 sources in the Fermi LAT First Source 
Catalog is made available by the Fermi Science Support Center at 
[http://fermi.gsfc.nasa.gov/ssc/data/access/lat/lyr_catalog/ 
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Figure 1. Smoothed Fermi-LKT count maps of front converted events in the range 2-200 GeV (units are 10"* counts deg"^), of the regions surrounding the 
remnants studied. The pixel binning is 0°.01, and the maps are smoothed with Gaussians of width 0°.2. In (a) Chandra X-ray emission map contours are overlaid. 
Black crosses indicate the direction of the detected OH (1720 MHz) maser emission. For CTB 37A, in (b), radio emission, as observed bv Whiteoak & Greet] 
(1996) in the Molonglo plane survey at 843 MHz, is overlaid as green contours. White crosses indicate the direction of the detected OH (1720 MHz) maser 
emission at velocities ~ —65 km s"' associated with CTB 37A, and the red crosses show the positions of maser emission at ~ -22 km s"', associated with 
G348.5-H0.0. The blue diamond indicates the position of CXOU J171419. 8-383023, and the cyan circle shows the Icr containment radius of HESS J1714-385. In 
(c), the green contours show the 20 cm emission toward 3C 391, as seen by the VLA FIRST Galactic Plane Survey, see White et al. (2005). Additionally, the cyan 
contours show the Chandra X-ray emission map. Yellow crosses indicate the direction of the two detected OH (1720 MHz) maser emission spots associated with 
3C 391. In (d), the green contours show the diffuse 20 cm emission toward G8. 7-0.1, as seen by the VLA FIRST Galactic Plane Survey. The black cross indicates 
the direction of the detected OH (1720 MHz) maser emission. The cyan ellipse contains the area where soft X-ray emission (0.1-2.4 keV) is observed with ROSAT, 
and the red circle indicates the approximate extension of HESS J1804-216.The black circle contains the radio extension of SNR G8. 31-0. 09, recently detected by 
la-ogan et al. (2006). The blue crosses indicate the position of Chandra sources CXOU J180432.4-214009 and CXOU J180441.9-214224 (Bamba et al. 200l 
iK^ galtsev et al. 2007a) . The magenta cross indicates the position of the young Vela-like pulsar B1800-21, found bv .Clifton & Lvne, ( ,1986,) ; ,Kargaltsev et all 
<2007bD . 



observations. iLazendic et all (120051) find that the overall X-ray 
emission of the remnant is best- fit by a nonequilibrium ioniza- 
tion plasma model with enhanced Si abundance, and temper- 
ature kT keV. The fit is significantly improved with the 
addition of a softer component of equilibrium thermal plasma 
with solar abundances and temperature kT w 0.8 keV. The Si 
overabundance suggests that the harder spectral component 
is emission from the supernova ejecta, and the softer compo- 
nent, with solar abundances, is related to shocked circumstel- 
lar material. From the volume emission measure these au- 
thors estimate the hydrogen density for the soft component to 

be «5 > cm"^, where /y is the volume filling factor, 

and assuming that the distance is 22 kpc and iie ~ I. Inn. The 
Chandra observations also revealed X-ray emission from a 
point source, CXOU J17 1801. 0-3726 17, inside the SNR shell, 
with X-ray characteristics similar to the compact central ob- 
ject (CCO) class, such as the source in Cas A, and unlike those 



for energetic pulsars (iPavlov et al.ll2004l) . 

The Fermi LAT count map around G349.7-I-0.2 is shown 
in Figure [ih, overlaid with the contours of the Chandra X- 
ray map. A bright unresolved 7-ray source is coincident with 
the position of the remnant, and with that of the source des- 
ignated as IFGL J1717. 9-3729 in the one-year Fermi LAT 
First Source Catalog. The significance of the detection, ob- 
tained from the evaluation of the test statistic in the region, 
is ^ 10(7. The residual test statistic map shows no evidence 
that the source is spatially extended. The spectrum of the 7- 
ray emission is shown in Figure |2^, and its shape suggests 
that the spectral index steepens above a few GeV. The best-fit 
model for the data yields a spectral index ofr = -2.10±0.11. 
The addition of an exponential cutoff is found to marginally 
improve the fit, the cutoff energy is estimated to be ^ 16.5 
GeV, and the spectral index is F = -1 .74 ± 0.37. Table[T|sum- 
marizes the characteristics of the best fit spectral model of the 
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Figure 2. Femii-LAT SED of the sources coincident with SNRs (a) G349.7+0.2, (b) CTB 37A, (c) 3C 391, and (d) G8.7-0. 1 . Statistical uncertainties are shown 
as black en'or bars, and systematic en'ors are indicated by red bars. The best-fit power-law models are shown (solid black hnes), as well as the spectra from 
TT^-decay models for proton distributions with energy cutoff at 100 TeV (dashed curves). For G349.7-I-0.2 and CTB 37A the spectra resulting from lower proton 
energy cutoffs, 160 Ge'V and 80 Ge'V respectively, are also shown (dotted curves). 



emission from this source, as well as its position. The lumi- 
nosity in the 0.1-100 GeV band of this source is estimated to 
be - 1.3x10^^- 



erg s 



, at a presumed distance of 22 kpc. 
3.2. CTB37A 



The SNR CTB 37 A (G348.5H-0.1) is located in a com- 
plicated region composed of three different remnants (Fig- 
ure [TJ)). Its radio emis sion overlaps that o f a second SNR, 
G348. 5-0.0 on the east (iKassim et al.lll991h . and a third rem- 
nant, CTB 37B, lies approximately 20' to the north. The ra- 
dio morphology of CTB 37 A appears as a well-defined partial 
shell towards the north and east, and an exte nded outbreak 
to the south. H 1 absorption measurements by ICaswell et al.l 
(119751) . toward both CTB 37A and CTB 37B, place con- 
straints on the distance to the CTB 37A/B complex, which 
is estimated to lie between 6.6 kpc and 11. 4 kpc, using a dis- 
tance to the Galactic center of 8.5 kpc. 

Various OH masers have been detected tow ard the interior 
of CTB 37A at 1720 MHz (iFrail et al.lll996h . Eight of these 
have velocities close to the upper velocity limit for this object 
given by the n eutral hydrogen measuremen ts (-65 km s~'). 
At this velocity. [Reynoso & MangumI (120001) observed, in the 
CO (7=1^0) transition at 1 15 GHz, three molecular clouds 
coincident with the maser spots and estimated the distance to 



th e clouds to be 11 . 3 kpc. 

lAharonian et al.l ( 120081) detected extended thermal X-ray 
emission from inside the NE shell of the remnant, using ob- 
servations with Chandra and XMM-Newton. This study found 
consistent fit results from the detectors on both X-ray obser- 
vatories, with plasma temperature T~ 0.8 keV and absorption 
A^H 3 X 10^^ cm~^. Additionally, these authors observed 
an extended non-thermal X-ray source, CXOU J171419.8- 
383023, toward the NW of the CTB 37 A complex (see Figure 
[T]3), which they suggest might be associated with the emis- 
sion of a pulsar wind nebula (PWN). A VHE (E >100 GeV) 
7-ray source, HESS J1714-385, is coincident with CTB 37A 
(Figure [TJ)), and it has been suggested this is a signature of 
the interaction between the remnan t's shock wave and the sur- 
rounding clouds of dense material (lAharonian et al.ll2008h . 

In Figure [T}^, the smoothed GeV map of the region around 
CTB 37A is shown. The 7-ray emission in the Fermi LAT 
range appears as an unresolved source whose centroid lies 
within the eastern radio shell of CTB 37A, with detection sig- 
nificance 17(7. This source is coincident with IFGL J1714.5- 
3830. The residual test statistic map gives no indication that 
this source is extended beyond the PSF of the instrument. The 
spectrum of this source is displayed in Figure together 
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Table 1 

Summary of spatial and spectral fit parameters from the Fermi-hAI data 







Spatial 




Spectral Fit 








R.A. 


Decl. 




F(O.l-lOOGeV) 








Name 


(deg) 


(deg) 


(deg) 


(10"' photons cm"^ s"') 


r 




TS Value 


G349.7+0.2 


259.47 


-37.50 


0.03 


0.58±0.11 


2.10±0.11 


1.12(5) 


108 


CTB 37A 


258.65 


38.52 


0.04 


1.36±0.15 


2. 19 ±0.07 


2.64(6) 


294 


3C 391 


282.26 


-0.92 


0.03 


1.58±0.26 


2.33±0.11 


0.80(6) 


140 


G8.7-0.1 


271.33 


-21.64 


0.03 


3.88 ±0.42 


2.40 ±0.07 


1.22(6) 


486 



with the best-fit models. The simple power-law model for the 
data yields a spectral index F = -2.19 ± 0.07. However, the 
fit is slightly improved by including an exponential cutoff at 
ficut ^ 4.2 GeV, with spectral index T = -1.46 ±0.32. Table 
[Usummarizes the result of the analysis of this source. The lu- 
minosity in the 0.1-100 GeV band of this source is estimated 
to be ~ 1 .3 X lO^'^erg s"' , at a presumed distance of 11 .3 kpc. 

3.3. 3C391 

3C391 (G31.9H-0.0) is a radio bright SNR with an irreg- 
ular spatial structure. Observations with the Very Large 
Array (VLAj^ in the radio band uncovered an elongated 
morphology, which exten ds from northwest to southeast 
jRevnolds & Moffet3ll993h . The shell surrounding the radio 
emission is enhanced toward the NW, and does not extend to 
the SE, where it appears that the SNR has broken out into a 
lo wer density regio n . 

iFrail et all (fT99^ observed two OH maser spots coincident 
with 3C 391 at 1720 MHz, with velocities (105 and 110 km 
s~') that a gree well with the neutr al hydrogen absorption ve- 
locity (Ra dhakrishnan et"aLlll972l) . The Hi absorption mea- 
surements indicate a distance of at least 7.2 kpc, and the ab- 
sence of absorption at negative velocities suggests an upper 
limit of 1 1 .4 kpc, using a distance to the center of the Galaxy 
of 8.5 kpc. The case for the interaction of 3C 391 with a 
molecular cloud, as suggested by the maser emission in the 
direction of the SNR, is further strengthened by the observa- 
tions of CO (y = 1 — > 0) emission in the region jWilner et alj 
[T998h . The CO maps show evidence that 3C 391 resides at 
the edge of a molecular cloud located to the NW of the SNR, 
and that the "breakout" morphology is produced as the SNR 
blast wave reaches lower de nsity regions oiitside the cloud. 

O bservations with A SCA jChen & Slanell2001l) and Chan- 
dra dChen et al.l 120041) show that the X-ray emission associ- 
ated with this object is thermal and clumpy in nature, and its 
brightness peaks well inside the radio shell. The plasma tem- 
perature derived from the X-ray observations is ~ 0.5 - 0.6 
keV, and the absorption column density is found to increase 
from SE to NW (from 0.6 to 4.1 x 10-^cm"-), which is con- 
sistent with the ex istence of molecular cloud to the northwest 
jChen et al] 120041) . These authors also discovered an unre- 
solved X-ray source on the NW border, with spectral charac- 
teristics similar to those of CCO class of objects. 

The smoothed photon count map from the ferm/-LAT data 
is shown in Figure [l];, and presents a strong case for an as- 
sociation between a GeV source and 3C 391. The unresolved 
7-ray source is detected with a 13cr significance and the peak 

' The VLA is operated by the National Radio Astronomy Observatory, 
which is a facility of the National Science Foundation, operated under coop- 
erative agreement by Associated Universities, Inc. 



of the test statistic map is 4' from the northwestern radio shell 
and the position of the closest maser spot. This source is 
coincident with IFGL J1849.0-0055. The SED of the emis- 
sion from this source is best described by a power law model, 
with spectral index T = -2.33 ± 0.1 1. The SED is shown in 
Figure |2]:, and spectral characteristics are presented in Table 
[U The luminosity of this source in the 0.1-100 GeV band is 
~ 0.9 X lO^^erg s"', at a presumed distance of 8 kpc. 

3.4. G8.7-0.1 (W30) 

The rather large SNR G8.7-0. 1, together with nine discrete 
Hn regions, forms the W30 complex. G8. 7-0.1 was identified 
as a Galactic SNR on the basis of the observation of a 45^ di- 
ameter area of extended non-t hermal radio emission jOdegardl 
[T986tlK assim & Weileri[T990l). 

iHewi tt & Yus ef-Zadehl (l2009h detected a single OH maser 
at 36 km s~', on the eastern edge of G8. 7-0.1. This relatively 
bright emission at 1720 MHz does not coincide with any of 
the compact radio sources, and hence they associate it with 
the SNR, and using the rotation curve of lFich et all (119891) 
they estimate a distance to the remnant of 4.5 kpc. An ex- 
tended 7-ray source, HESS J1804-216, has been detected to 
the southwest of the remnant (Figure[T]i), which could be fur- 
ther e vidence of shock intera ction with clouds of dense ma- 
terial jAharonian et alj 120061) . Several other sources in the 
field have been suggested as possible associations with the 
source HESS J 1804- 2 16, and they are shown in Figure [TB 
jCHfton & Lvne|1986l:lBro gan et al. 2006; Bamba et al. 200l 
iK argaltsev et al. 2007a b; Higa shi et al.l l2008). 

iFinlev & Ogelman (1994) detected diffuse thermal X-ray 
emission from the northern half of the remnant in the 0.1-2.4 
keV band with ROSAT, and using Sedov analysis estimated 
the distance to the remnant to be 3.2 kpc < d < 4.3 kpc. 
The estimated plasma temperature from these X-ray obser- 
vations is r « 0.3-0.7 keV, the absorption column density is 
A^H ^ 1 .2- 1 .4 X 10^^ cm~^, and the derived ambient density is 
«o ^ 0.02-0.04 cm"-'. These authors also determined the dis- 
tance to the Hll regions associated with re mnant to be ~ 4.8 
kpc, u sing the kinematic measurements of iKassim & WeileJ 
(119901) and more recent Galactic rotation models. 

The 7-ray photon count map in the 2-200 GeV range of 
the crowded vicinity of the W30 complex is shown in Fig- 
ure[T]i. The 7-ray emission in the region appears as a source 
centered on a position inside the radio boundary of G8.7-0. 1 . 
This detection is coincident with IFGL J1805. 2-2137. The 
residual test statistic map does not rule out the possibility 
that this source is extended. The significance of the detec- 
tion is ^ 22(T. The spectrum of this source, shown in Fig- 
ure |2}l, seems to steepen above 2 GeV, and flatten at > 10 
GeV. For a single power-law model the best-fit spectral index 



6 



Castro & Slane 



is r = -2.40 ± 0.07. The spectral characteristics of this region 
are summarized in Table [T|. The luminosity in the 0.1-100 
GeV band of this source is 
distance of 4.5 kpc. 



1.1 X 10 erg s , at a presumed 



Table 2 

Model parameters, distances, and density estimates 



4. DISCUSSION 

The detection with the Fermi-LAT of unresolved 7-ray 
sources in the direction of G349.7+0.2, CTB 37A, and 3C 
391, as well as a possibly extended region of emission co- 
incident with G8. 7-0.1 is reported in this work. The best-fit 
positions and spectral characteristics of these sources are sum- 
marized in Table [T] 

The shapes of the spectra of the sources coincident with 
G349.7-0.5, 3C 391 and G8.7-0. 1 do not seem con sistent with 
what has been observed from pulsars in 7-rays. lAbdo et aP 
(|2009b) found the SED of pulsars in the Fermi LAT band to be 
best described by power-law distributions, with exponential 
cutoffs in the energy range 1-6 GeV (the only exception, PSR 
J 1833- 1034 with cutoff energy E^ut ^ 8 GeV, was reported 
to be poorly measured). The spectral characteristics of the 
sources coincident with 3C 391 and G8. 7-0.1 do not suggest 
an exponential cutoff at all, and the spectrum of G349.7-I-0.2 
is very modestly improved by a cutoff at E^ut = 16.5 ±1.9 
GeV, which is very high for a 7-ray pulsar. Additionally, 
there are no pulsars from the Australian Telescope National 
Facility (ATNF) Pulsar C atalogue^ coincident with this source 
(iManchester et al.ll2005h . 

The spectrum of the source coincident with CTB 37A is 
well fit by a power law with an exponential cutoff at energy 
Ecut ^ 4.2 GeV, and hence does not rule out a pulsar hypothe- 
sis. However, there are no ATNF pulsars within 20' radius of 
this source. 

Given the lack of evidence for significant contribution by 
a pulsar, and the known presence of maser emission for each 
remnant, an appropriate hypothesis appears to be that the 7- 
ray sources are products of SNR interactions with molecular 
clouds. To investigate properties of the expected spectra, we 
have considered the scenario in which the 7-rays are produced 
by 7r° decay of accelerated hadrons. In order to reproduce the 
observed spectra, w e have used a model based on that from 
iKamae et al.1 (l2006l) using a scal ing factor of 1.85 for helium 
and heavier nuclei ( iMoril |2009|) . and adopted a power law 
model with an exponential cutoff for the energy distribution 
of the protons: 



dNp 
dE„ 



1 TeV 



exp 



(1) 



where £0 is the proton energy cutoff. Table |2] presents sets 
of representative model parameters that adequately reproduce 
the observed 7-ray spectra (see Figure 2). While the spectra 
of SNRs 3C 391 and G8.7-0.1 are well matched by models 
with proton distributions with spectral indices F w 2.4 and 
high-energy cutoff, £0 > 100 TeV, the spectra of G349.7H-0.2 
and CTB 37 A require harder proton distributions, F « 1.7, 
and £() ~ 10~ GeV. The model provides an estimate of the 
density of the material the SNR shock wave is impacting, 
given reasonable assumptions about the shock compression 
ratio, r, and the fraction e of the total supernova explosion 
energy E that is converted into cosmic ray energy. Here we 

* The positions and characteristics of the 1509 pulsars in ATNF 
catalogue, as well as further documentation, can be accessed at 
[http://www.atnf.csiro.au/re search/pulsar/psrcat 





Model Parameters 






Object 


r„ 


Eo 




"x 


d 






(TeV) 


(cm-3) 


(cm-3) 


(kpc) 


G349.7+0.2 


1.7 


0.16 


65 


2.5 


22 


CTB 37A 


1.7 


0.08 


37 




11.3 


3C 391 


2.4 


100 


28 


1 


g 


G8.7-0.1 


2.45 


100 


18 


0.03 


4.5 



"Densities derived using the 7-ray spectrum obtained from the analy- 
sis of Fermi-LNI observations. 

''Densities derived from X-ray observations (see text for references). 
No estimate of the ambient density of SNR CTB 37A is provided in 
Aharonian et al. (2008). 



assume r = 4, e = 0.4, and E = 10^' erg; we use distances de- 
rived from observations at other wavelengths for each SNR. 
In Table |2l estimates of the densities of the 7-ray emitting 
material are compared to the densities derived fro m X-ray ob- 
serva t ions in other studie s (Lazendicet al. 2005; Chen et all 
l200llFiniev & OgelmanlfT994 . 

The values obtained from the Tr^-decay model fits to the 
7-ray spectra of the SNRs are much higher than those from 
X-ray observations. It is possible that, upon interacting with 
dense ambient material, instabilities in the postshock flow re- 
sult in considerable clumping of material. The denser regions 
would then not be expected to yield significant X-ray emis- 
sion, although radiative shocks would then be expected to 
produce optical emission from these regions. If the filling 
fraction of the clumped material is large, then the densities 
inferred for the X-ray emitting material are underestimates, 
potentially leading to more reasonable values for the density 
contrast between clump and interclump regions. 

Alternatively, the energetic particles that escape the particle 
acceleration region at the SNR shock could stream ahead of 
the shock and encounter clouds of molecular material, where 
the much higher density could dramatic a fly enhance the 7-ray 
emission (lAharonian & AtoyanI Il996t iGabici & Aharonianl 
l2007h . This scenario has been proposed as a possible ex- 
planation of the observed 7-ray fluxes in the di rection of 
some Galactic S NRs, including RX J 17 13.7-3946 (iButt et al.l 
l200lb and W28 dFuiita et al.ll2009l) . [Lee et al.l (|2008|) model 
the interaction of particles accelerated by the diffusive shock 
of an SNR with arbitrary matter distributions outside the 
remnant blast wave, including an appropriate treatment of 
the escape and diffusion of th ese particles. These au- 
thors, as well as other works dOabici & Aharonianl [20071: 
[Rodriguez Marrero et al.|[2008l: [Gabici et al.[[2008[) ," show that 
clouds of dense material in close proximity to SNRs can have 
an important, and even dominant contribution to the 7-ray 
flux from these regions. It is important to note, however, that 
the escaping particles are predominantly from the high-energy 
portion of the spectrum; the observed spectra, which extend 
to low energies, may be problematic for such a scenario. 

An alternative scenario is that the 7-ray emission coinci- 
dent with these sources is of leptonic origin. Bremsstrahlung 
and inverse Compton (IC) radiation from the nonthermal elec- 
tron population have been pr oposed as producti on mecha- 
nisms of 7-rays in SNRs (e.g. [Bvkov et al.|[2000l) , but these 
models present certain difficulties as well. Non-thermal 
bremsstrahlung can be the dominant emission process in this 
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band for electron-to -proton ratios, ae/up, greater than ^ 0.2 
jGaisser et al.l 119981) . However, cosmic -ray abundance ra- 
tios suggest flg/flp ~ 10~^, and models for 7-ray emission 
from other SN Rs (e.g RX J1713. 7-3946) favor even sm aller 
values (Morlin o et al.ll2009t IZirakashvili & Aharoniarill2010l: 
lElhson et al.. .20101) . making the bremsstrahlung-dominated 
model appear less likely. The IC-dominated scenario, on 
the other hand, is difficult to reconcile with the observed 
characteristics of the sources studied in this work. The ex- 
pected spectral index for IC emission (F « 1 .5) is consider- 
ably harder than the values obtained here (F > 1 .7). 

In summary, using observations with the FermZ-LAT, we 
have revealed 7-ray emission from four SNRs known to be 
interacting with molecular clouds. In each case the inferred 
density is high, assuming a hadronic origin for the 7-rays, 
consistent with the presence of dense material indicated by 
maser emission. However, estimated densities of the ambi- 
ent material derived under simple geometric assumptions are 
in considerable excess of those derived from the hot X-ray 
gas observed in the SNRs. More detailed modeling of the X- 
ray and 7-ray emission under a variety of scenarios for the 
ambient medium is required to address the emission mecha- 
nism more f ully. These results ar e consistent with the scenario 
presented bv lAbdo et all (I2009dl) for SNR W5 IC, and further 
studies of other such systems and their GeV-scale 7-ray emis- 
sion with the Fermi-hP^ are of particular interest. 

The authors thank Yasunobu Uchiyama and Stefan Funk 
for their careful review of this work, and Don Ellison, Keith 
Bechtol, Joshua Lande, Andrei Bykov, Luke Drury, and Felix 
Aharonian for helpful discussions. This work was partially 
funded by NASA contract NAS 8-03060 and Fermi grant 
NNX09AT68G. RO.S. acknowledges the KITP in Santa Bar- 
bara, where work on this project was begun while participat- 
ing in a KITP program. 
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